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We present an exact method based on a hyperspherical harmonic expansion to study 
systems made of quarks and antiquarks of the same flavor. Our formalism reproduces 
and improves the results obtained with variational approaches. This analysis shows that 
identical-flavor four-quark systems with non-exotic 2 ++ quantum numbers may be bound 
independently of the quark mass. + ~ and l + ~ states become attractive only for larger 
quarks masses. 

1. Introduction 

The understanding of few-body systems relies in our capability to design methods for 
finding an exact or approximate solution of the iV— body problem. The solution of any 
few-particle system may be found by means of an expansion of the trial wave function in 
terms of hyperspherical harmonic (HH) functions. The idea is to generalize the simplicity 
of the spherical harmonic expansion for the angular functions of a single particle motion 
to a system of particles by introducing a global length p, called the hyperradius, and 
a set of angles, Q. For the HH expansion method to be practical, the evaluation of the 
potential energy matrix elements must be feasible. The main difficulty of this method is to 
construct HH functions of proper symmetry for a system of identical particles. This may 
be overcome by means of the HH formalism based on the symmetrization of the N— body 
wave function with respect to the symmetric group using the Barnea and Novoselsky 
algorithm 1 . 

The recent discoveries of several meson-like resonances whose properties do not fit 
into the predictions of the naive quark model, has reopened the interest on the possible 
role played by non-gg configurations in the meson spectra. Among them, the possible 
existence of bound four-quark states (two quarks and two antiquarks) has been suggested 
in the low-energy hadron spectroscopy. Four-quark bound states were already suggested 
theoretically thirty years ago, both in the light-quark sector by Jaffe[2] and in the heavy- 
quark sector by Iwasaki[3]. 

In this work we present a general study of four-quark systems of identical flavor in 
an exact way. For this purpose we have generalized the HH method, widely used in 
traditional nuclear physics for the study of few-body nuclei, to describe systems made 
of quarks and antiquarks. This generalization presents two main difficulties, first the 
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simultaneous treatment of particles and antiparticles, and second the additional color 
degree of freedom. 

2. General formulation of the problem 

The system of two quarks and two antiquarks with the same flavor can be regarded as 
a system of four identical particles. Each particle carries a SU (2) spin label and a SU(3) 
color label. Both quarks and antiquarks are spin 1/2 particles, but whereas a quark color 
state belongs to the SU(3) fundamental representation [3], an antiquark color state is a 
member of the fundamental representation [3]. The four-body wave function is a sum of 
outer products of color, spin and configuration terms 

|0) = | Color) | Spin) | R) (1) 

coupled to yield an antisymmetric wave function with a set of quantum numbers that 
reflects the symmetries of the system. These are the total angular momentum quantum 
number J, its projection J 2 , and the 577(3) color state G (labeled as G to avoid confusion 
with the charge conjugation quantum number), which by assumption must belong for 
physical states to the S77(3) color singlet representation. Since QCD preserves parity, 
parity is also a good quantum number. Another relevant quantum number to the system 
under consideration it is the C— parity, C, i.e., the symmetry under interchange of quarks 
and antiquarks. 

To obtain a solution of the four-body Schrodinger equation we eliminate the center of 
mass and use the relative, Jacobi, coordinates ffi, if 2 , . . . , fjA-i- Then we expand the spatial 
part of the wave-function using the HH basis. In this formalism the Jacobi coordinates 
are replaced by one radial coordinate, the hyperradius p, and a set of (3 A — 4) angular 
coordinates Qa- The HH basis functions are eigenf unctions of the hyperspherical part 
of the Laplace operator. An antisymmetric A-body basis functions with total angular 
momentum J a, J a , color Ga and C-parity C, are given by, 
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where (p\n) = R n (p) are the hyperradial basis functions, taken to be Laguerre functions. 

{Qa^aLaMaVaYa-^a) = y [ £ ] A L A M A T A Y A _ iaA ^A) (3) 

are HH functions with hyperspherical angular momentum K = Ka, and orbital angular 
momentum quantum numbers (La, Ma) that belong to well-defined irreducible representa- 
tions (irreps) Ti G T 2 . . . G of the permutation group-subgroup chain S\ C S2 ■ ■ ■ C Sa, 
denoted by the Yamanouchi symbol [r^, Y A _i] = \T A , r^_ 1? . . . , Ti\. The dimension of the 
irrep T m is denoted by |T m | and A rA y A1 is a phase factor. Similarly, the functions 

{s' 1 .. 8 ' A ,g 1 ..g A \S A S' A G A Cf A ,Y A . 1 p A ) = x^^y^Jsl.s^, 9l ..g A ) (4) 

are the symmetrized color-spin basis functions, given in terms of the spin projections 
(s 2 ) and color states (&) of the particles. The quantum numbers ot A , and (3 A are used to 
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remove the degeneracy of the HH and color-spin states, respectively. For the construction 
of the symmetrized HH basis we use the algorithm of Barnea and Novoselsky PQ, which 
utilizes the group of kinematic rotations. For the color-spin subspace, we use a method 
to transform the standard basis into a symmetrized color-spin basis with well defined 
color and C— parity. The technical steps of such construction are detailed in Ref. [3]. 
The calculation of the Hamiltonian matrix-elements between the antisymmetric basis 
functions, Eq. (j2J), is almost the same as in nuclear physics, replacing isospin by color 

Table 1 

cccc masses for the maximum value of K used, E(K max ), and using the extrapolation, 
E(K = oo), compared to the corresponding threshold, M±M 2 and T(M l5 M 2 ). The 
subindex stands for final state relative angular momentum. The value of A for each 
state is also given. Energies are in MeV. 
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3. Results 

We have applied our method to study cccc states. The calculation has been done up to 
the maximum value of K within our computational capabilities (K max ). To analyze the 
stability of these systems against dissociation through strong decay, parity (P), C— parity 
(C), and total angular momentum (J) must be conserved. 

We show in Table ^ all possible J PC quantum numbers with L = 0. We also indicate 
the lowest two-meson threshold for each set of quantum numbers. Let us note that 
the convergence of the expansion in terms of hyperspherical harmonics is slow, and the 
effective potential techniques [5 J do not improve it. To obtain a more adequate value for 
the energy we have extrapolated it according to E(K) = E{K = oo) + a/K b , where 
E{K = oo), a and b are fitted parameters. The values obtained for E(K = oo) are stable 
within ±10 MeV for each quantum number. Four-quark states will be stable under strong 
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interaction if their total energy lies below all possible, and allowed, two-meson thresholds. 
It is useful to define A = E(K = oo) — T(M 1; M 2 ) in such a way that if A > the four- 
quark system will fall apart into two mesons, while A < will indicate that such strong 
decay is forbidden and therefore the decay, if allowed, must be weak or electromagnetic, 
being its width much narrower. 

A first glance to these results indicates that only three sets of quantum numbers have 
some probability of being observed. These are the J PC = + ~, and J PC = 2 ++ , 
which are very close to the corresponding threshold. It is interesting to observe that 
the quantum numbers + ~ correspond to an exotic state, those whose quantum numbers 
cannot be obtained from a qq configuration. 

To analyze whether the existence of bound states with exotic quantum numbers could 
be a characteristic feature of the heavy quark sector or it is also present in the light sector 
we have calculated the value of A for different quark masses. We have obtained that 
only one of the non-exotic states, the 2 ++ , becomes more bound when the quark mass is 
decreased, A pa —80 MeV for the light quark mass. The l +_ and ++ states, that were 
slightly above threshold in the charm sector, increase their attraction when the quark 
mass is increased and only for masses close to the bottom quark mass may be bound. 
With respect to the exotic quantum numbers, the negative parity and 1~ + are not 
bound for any value of the quark mass. Only the + ~ four-quark state becomes more 
deeply bound when the constituent quark mass increases, and therefore only one possible 
narrow state with exotic quantum numbers may appear in the heavy-quark sector. This 
state presents an open P— wave threshold only for quark masses below 3 GeV. 

Let us note that since the heavy quarks are isoscalar states, the flavor wave function of 
the four heavy-quark states will be completely symmetric with total isospin equal to zero. 
Therefore, one should compare the results obtained in the light-quark case with a com- 
pletely symmetric flavor wave function, i.e., the isotensor states. There are experimental 
evidences for three states with exotic quantum numbers in the light-quark sector. Two 
of them are isovectors with quantum numbers J PC = 1~ + named 7Ti(1400) and 7Ti(1600), 
and one isotensor J PC = 2 ++ , the X(1600). Taking the experimental mass for the thresh- 
old T(Mi, M 2 ) in the light-quark case together with the values obtained for A, one can 
estimate the energy of these states, being M{2 ++ ) pa 1500 MeV and M(l~ + ) pa 2900 
MeV. The large mass obtained for the 1~ + four-quark state makes doubtful the identifi- 
cation of the 7r!(1400) or the 7Ti(1600) with a pure multiquark state, although a complete 
calculation is needed before drawing any definitive conclusion. Concerning the V(1600), 
being its experimental mass 1600±100 MeV, a tetraquark configuration seems likely. 
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